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Abstract 

We examine the effects of spatial inhomogeneities on irrotational anisotropic cosmolo- 
gies by looking at the average properties of anisotropic pressure-free models. Adopting the 
Buchert scheme, we recast the averaged scalar equations in Bianchi-type form and close the 
standard system by introducing a propagation formula for the average shear magnitude. We 
then investigate the evolution of anisotropic average vacuum models and those filled with 
pressureless matter. In the latter case we show that the backreaction effects can modify 
the familiar Kasner-like singularity and potentially remove Mixmaster-type oscillations. The 
presence of nonzero average shear in our equations also allows us to examine the constraints 
that a phase of backreaction-driven accelerated expansion might put on the anisotropy of 
the averaged domain. We close by assessing the status of these and other attempts to define 
and calculate 'average' spacetime behaviour in general relativity. 

PACS Numbers: 04.20.-q, 04.40.-b, 98.80.-k, 98.80.Jk 

1 Introduction 

The question of whether an inhomogeneous model of the universe behaves on the average like a 
homogeneous solution of Einstein's equations is a long-standing problem, directly related to the 
question of how to average general-relativistic spacetimes [l]-[3]. The problem lies in the generic 
non-linearity of general relativity and in the difficulty of establishing a unique and unambigu- 
ous way of averaging the spacetime metric without throwing away crucial information in the 
averaging process. A number of cosmologists have argued that the averaging problem may be 
central to our understanding of the recent expansion history of the visible universe. When it was 
suggested, recently, that structure formation could be responsible for the observed acceleration 
of the universe [4] -[7], averaging techniques were employed to study the dark energy problem in 
terms of kinematic backreaction effects from spatial inhomogeneities [8, 9]. Although the idea 
of backreaction has been criticized, primarily on the basis of current observations [10]- [13], it is 
attractive in principle because it can solve the dark energy and the coincidence problems without 
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appealing to a cosmological constant, speculative quintessence field, or non-linear corrections in 
Einstein's theory of gravity. 

General relativity has been applied to a range of averaging scales depending on the nature 
of the physical system under study. Cosmology applies to the largest astronomical scales, which 
typically extend over a significant fraction of the Hubble radius. Nevertheless, the major issue 
of defining a suitable averaging process remains open. When employing averaging techniques 
to study the backreaction of cosmological perturbations, one usually averages scalars like the 
density or the volume expansion rate. These become the scale-dependent parameters that are 
compared with observations. Averaging such scalars over a Riemannian spacelike region gives 
the effective dynamical sources that an observer would expect to measure. Thus, the averaged 
equations isolate an explicit source term that is commonly referred to as the 'backreaction' term. 
The latter quantifies the deviation of the average model from a given 'fitting' model, which is 
usually the Friedmann-Robertson- Walker (FRW) universe, though not always [14]. 

Here, we do not use a Friedmannian universe, but an anisotropic Bianchi-type cosmology that 
contains the FRW models as subcases. We consider vacuum, or dust-dominated, irrotational 
Bianchi models with zero or isotropic local spatial curvature, aiming to describe the effective 
evolution of their averaged counterparts. We do this by writing the standard average formulae in 
Bianchi form (i.e. with explicit shear terms in the left-hand side) and by closing the system with 
a propagation equation for the mean shear magnitude. Within this Bianchi-type interpretation 
of the averaged equations the backreaction term is expressed as a difference between the mean- 
square fluctuations in the expansion and shear scalars. For a vacuum background with Euclidean 

3- spaces, we find that the source term vanishes and the average equations retain their local form. 
When pressureless matter is introduced, however, the backreaction effects can dominate the 
early evolution of averaged Bianchi / models and prevent the familiar Kasner-type singularity 
from forming. Assuming axial symmetry, we examine the 'bounce conditions' that avoid the 
initial singularity for different types and degrees of anisotropy. We also consider the accelerated 
expansion of these models and how this might constrain their average shear anisotropy. 

2 The averaging rule 

The literature contains more than one averaging method [16]-[22] that has been applied to 
cosmology (see [23] and more recently [8, 9] and [24]- [27]). Here, we will follow the Buchert 
approach which builds on the Newtonian theory, where spatial averaging is relatively simple, 
before extending to general relativity by confining to scalar variables [20]-[21]. Consider a gen- 
eral spacetime filled with irrotational pressureless matter moving along geodesic worldlines with 

4- velocity u a (i.e. u a u a = — 1 and u a = u b \/i,u a = 0). Given a foliation of flow-orthogonal hyper- 
surfaces, the metric of the 3-dimensional space is = g a ^ + u a Uf,. Then, the spatial averaging 
of an arbitrary scalar field 4> = 4>(x a ) over a simply-connected domain V is a covariant operation 
defined by the averaging operation {4>)t> = Jj? fiJd^x /Vr>, where the angled brackets indicate 
spatial averages, Vp is the volume of the domain in question, and J = \J det(/i a ;,) [20]. We also 
note that Greek indices take the values 1, 2, 3 and Latin ones run from to 3. Following [8], it 
helps to define an effective expansion scale factor by ax> = (Vp/V^) 1 / 3 , with Vxi i representing 
the volume of the initial domain. Since the volume Vx> is generally time dependent, the average 
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of the time derivative of a locally-defined scalar differs from the time derivative of the average. 
This non-commutativity between spatial averaging and temporal evolution can be formulated 
in a simple mathematical rule. For a rest-mass preserving domain, the commutation between 
the volume-averaging and the time-evolution operators leads to [20, 21] 

Wp-Wd = W>d-(«d, (1) 

where is the local volume expansion rate between neighbouring worldlines. This scalar coin- 
cides with the trace of the local expansion tensor and also with that of the extrinsic curvature 
tensor. The right-hand side of the above rule generally introduces source terms to the averaged 
equations which can be interpreted as backreaction effects due to the averaging process. Such 
kinematic backreaction has recently been considered as a possible solution to the dark energy 
problem. We note that the appearance and the nature of the aforementioned source terms de- 
pends on the form of the averaged equations. Generally, the less information these equations 
contain the simpler (if any) the source terms are. 

3 Averaged anisotropic cosmologies 
3.1 The local evolution 

Consider an irrotational, anisotropic, spacetime and introduce a family of (fundamental) ob- 
servers living along worldlines tangent to the future-directed timelike eigenvectors of the Ricci 
tensor. This achieves a unique local threading of the spacetime into time and 3-dimensional 
space (e.g. see [28]). For a pressure-free medium, the local shear tensor propagates as [29, 30] 

<j a b = -®v a b ~ K(ab) , (2) 

where is the traceless part of the locally-defined 3-Ricci tensor. The latter describes the 

geometry of the spatial hypersurfaces orthogonal to u a and satisfies the Gauss-Codacci equation. 
Hereafter, we will only consider spacetimes that obey the constraint lZ a ba ab = 0. For example, 
this condition holds when the curvature of the spatial sections vanishes, or when it is isotropic. 
Then, the local expansion is monitored by the following three scalar propagation equations 

p + 6p = 0, <j + Oa = and + 3-G 2 + \p + 2a 2 - A = , (3) 

supplemented by the constraint 

\K- p+\<d 2 - a 2 - A = 0. (4) 

Here, p is the matter density, G is the expansion scalar, a is the shear magnitude, TZ is the Ricci 
scalar of the spatial sections, and A is the cosmological constant. 

We remind the reader that the adopted averaging scheme requires zero rotation to guarantee 
the existence of integrable spatial hypersurfaces. Also, the method has been applied to study 
the backreaction effects in pressure-free cosmological environments. Here, in addition to zero 
rotation and pressure, we assume that lZ a i,a ab = 0. For all practical purposes the latter is a very 
mild constraint and still covers all (non-rotating and pressure-free) spacetimes with flat spatial 
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sections or isotropic 3-curvature. These include those Bianchi models that contain the FRW 
universes as subcases. Primarily, condition 1Z a t,cr ab = allows us to reduce Einstein's equations 
to a closed set of four scalar formulae. Previous work on the backreaction effects of spatial 
averaging did not need a shear propagation equation and the associated sets of scalar equations 
were not closed (in reminiscent of the Newtonian 'gravitational paradox' - see [31]). 

3.2 The average evolution 

Because spatial averaging and time evolution do not commute, the averaged equations contain 
extra terms, compared to their local counterparts. By giving to the average formulae the form 
of a chosen set of local equations, one can interpret these extra terms as the backreaction effect 
of the averaging process. Here, we will write the averaged formulae in Bianchi-type form by 
allowing for explicit shear terms in their left-hand side. Thus, applying law (1) to the matter 
density, the shear magnitude, and the volume expansion, we may use Eqs. (3a)-(3c) to obtain 

( P y + (e)(p) = o, ( a y + (e)(a)=o (5) 

and 

(©)' + U®) 2 + UP) + 2<^> 2 - A = !<(© - <©» 2 > - 2<(<7 - (a)) 2 ) , (6) 
respectively. Also, the volume- averaged counterpart of (4) reads 

- (p) + i<e> 2 - (a) 2 - A = 4<(6 - <9» 2 > + <(<x - (a)) 2 ) , (7) 

where the right-hand sides of (6) and (7) involve the identity ((4> — (<ft)) 2 } = (4> 2 ) — (4>) 2 ■ The 
above set is a system of four effective Einstein equations for spatially averaged scalar variables 
in inhomogeneous irrotational universes with zero fluid pressure and TZ ab a ab = 0. Within this 
environment Eqs. (5)-(7) are exact and there is no need to assume that the inhomogeneity and 
the anisotropy are small perturbations [8, 9]. 

Comparing Eqs. (5) to (3a) and (3b) we notice that, despite the non-commutativity between 
volume averaging and time evolution, the averaged formulae have preserved the form of their 
local counterparts. This means that 

(p) oc a^, 3 and {a) oc a^, 3 , (8) 

where a© is the average scale factor smoothed out over the domain V (with ax>/ax> = (0)/3). 
However, the averaging process has led to extra terms in the right-hand sides of (6) and (7), 
collectively given by the domain-dependent scalar 

Q = 2 <(0 _ (6)) 2 ) _ 2{{a _ {a)) 2 } = _2 ((@) 2 _ (@2)) + 2 ((a) 2 _ ^ (g) 

We interpret the above as the backreaction effect of spatial averaging on pressure-free Bianchi I 
models, since in its absence the averaged formulae recover the form of their local counterparts 
(compare to Eqs. (3), (4)). Note that Q can be either positive or negative, and its value is 
decided by the amount of kinematic anisotropy. In particular, Q is positive when the shear 
mean-square fluctuation is less than a third of that in the volume expansion. 
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3.3 The kinematic backreaction term 



The form of the Q-term depends on the specifics of the adopted cosmology. For instance, in 
empty models with zero cosmological constant and 3-curvature, the local shear and the expansion 
are related by 0/\/3 = a (see Eq. (4) earlier). This immediately ensures that (0)/\/3 = (cr), 
and subsequently that (0} 2 /3 = (cr) 2 . Hence, following (9), the backreaction term vanishes 
identically, and Eqs. (6), (7) maintain the form of their local counterparts. 

When there is matter or nonzero 3-curvature, the backreaction effects do not generally vanish. 
In those cases we can monitor the dynamical evolution of Q via a consistency /integrability 
condition. Taking the time derivative of (7) and then using (5), (6) we obtain 

Q + 2(@)Q = -(ny-K@)(n). (10) 

It follows that in the absence of pressure and rotation, and also for TZ a i ) a ab = in our case, the 
evolution of the backreaction term depends on that of the average 3-Ricci scalar. When 1Z = 
or (1Z) oc (a)~ 2 , for example, we find that 

Qota v 6 . (11) 

Overall, a positive Q mimics the dynamics of a dark energy component and tends to accelerate 
the expansion, whereas a negative one acts as an effective dark matter source (see Eq. (6)). In 
particular, a negative Q behaves like an effective fluid with a 'stiff' equation of state (e.g. a free 
scalar field) or like an additional shear source. 

4 Backreaction effects 

4.1 Backreaction and Kasner-type singularities 

The scalar Q carries the collective backreaction effects from spatial inhomogeneities. This could 
mean a considerable loss of information, relative to that contained in the local equations. For 
instance, the sign of the backreaction term is crucial but unknown. There is an advantage, 
however, because the compactness of Q makes it an easier to handle quantity. Thus, a posi- 
tive Q was recently suggested as potential cause of late-time acceleration in FRW-type average 
cosmologies. Here we will consider average Bianchi-type models and examine whether anal- 
ogous backreaction effects can modify (or even remove altogether) standard features of their 
conventional counterparts. 

Integrating (3a) and (3b), we find that p oc a~ 3 and a 2 oc a~ 6 , with the latter holding for 
TZ ab a ab = 0. This means that shear anisotropies dominate the early evolution of these models, 
which approach the Kasner-type singularities of the vacuum Bianchi / spacetimes (with A = 0). 
This singular behaviour is not guaranteed for the averaged counterparts of these cosmologies 
because of the backreaction term and particularly because of (11). Given that both (a) 2 and Q 
propagate as a^> 6 while (p) oc a^> 3 , the early stages of the averaged kinematics is dominated by 
the sum Q-2(a) 2 = 2((6 - (9)) 2 )/3 - 2(a 2 ). When the latter is negative and (6) < initially, 
the collapse will proceed unimpeded and (0) — ► -co within a finite amount of time. Otherwise 
the averaged model should bounce before reaching the singularity. 
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Additional information is obtained by expressing Q in terms of the averaged individual 
expansion rates. Assume, for example and also for the sake of simplicity, that the inhomogeneous 
spacetime is axially symmetric and has zero 3-curvature. Then 6 = 2a + (3, p = a(a + 2(3) and 
a = (a — (3)/y/3. In this notation a = a/a and [3 = bjb are the expansion rates along the two 
main axes, with a and b representing the associated scale factors of the axisymmetric Bianchi / 
universe [32, 33]. Within this environment, expression (9) reads 

Q = 2((a 2 )-(a) 2 )+A((aP)-(a)((3)) , (12) 

allowing for Q to be positive, negative, or zero. 1 Hence, if the expansion is nearly isotropic 
(i.e. for a ~ (3), the bounce condition Q > 2(a) 2 is satisfied when ((a — (a)) 2 } > 0. For 
highly anisotropic expansion, with a » /3, we have (a) 2 ~ (a) 2 /3. In this case, we will have a 
bounce when the mean-square fluctuation of a is greater than (a) 2 /3. Finally, for a <C (3, the 
backreaction effects can prevent the (0) — > — oo singularity provided that (a(3) — (a)((3) > ((3) 2 /6. 

It is well known that the Mixmaster behaviour of the Bianchi IX models consists of a 
sequence of chaotically alternating Kasner-like time intervals where the collapsing and expanding 
scale factors interchange (e.g. see [32] and also [34, 35]). The possibility that backreaction effects 
due to spatial averaging could alleviate the Kasner-type singularities of the Bianchi I models, 
raises the question as to whether the same effects could also modify the Mixmaster regime of 
the Bianchi IX spacetimes. The possible role of stochastic spatial fluctuations was discussed 
in [36]. It was argued there that the cumulative stochastic effects from nearby gravity- wave 
perturbations could transform the local Einstein equations into stochastic differential equations 
with significant damping effects on the chaotic oscillations. Here, we see the cumulative effect 
of the stochastic inhomogeneities mimics a scalar field of the sort that is known to stop the 
continuation of Mixmaster oscillations on approach to the initial singularity (if it exists) [37]. 
Based on these we conjecture that, within the framework of the Buchert averaging scheme, 
Mixmaster oscillations will be unstable at early times. 

4.2 Accelerated expansion and average anisotropy 

Let us suppose that the aforementioned backreaction effects lead to the acceleration of the 
averaged domain. The presence of nonzero mean shear in the equations allows us to consider 
the constraints imposed on the average anisotropy of the region. Using the scale-dependent 
average scale factor ax>, defined so that d^/a© = (B)/3, the averaged Raychaudhuri equation 
(see (6)) takes the form 

3^ = -l(p)-2(a) 2 + A+Q. (13) 

Here, the backreaction term can be interpreted as an effective energy-density term or an effective 
shear term. 2 In any case, and in the absence of a cosmological constant, the volume expansion 

1 When a = — 2/3 expression (12) ensures that Q = and therefore the absence of any averaging backreaction 
effects. This is expected since a = —2/3 means p = and we know that Q — in the vacuum models (see § 3.3). 

2 One can introduce a chosen 'background' to the average equations by defining a suitable reference model. 
Thus, the time-dependent expansion and shear fields, H = H(t) and E (, = E (,(t) respectively, define a homo- 
geneous background into the averaged formulae [20]. Inhomogeneous deviations from these reference fields are 
considered by setting O = H + 6 and a a b = S a 6 + a a b- This type of decomposition is not essential for our 
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of the averaged domain will accelerate provided that 

Q>±(p) + 2(v) 2 , (15) 

which here holds for zero or isotropic 3-curvature (strictly speaking as long as TZ a i,a ab = 0). 
We note that if Q is negative, namely if the mean square fluctuation in the shear is greater 
than a third of the mean square fluctuation in the volume expansion (see (9)), the domain will 
decelerate instead of accelerating. We have no information on the sign of the backreaction term 
however. Because of this ambiguity, early on, kinematic backreaction was also suggested as an 
effective dark-matter source [38]. 

Putting these matters aside, we will use our Bianchi-based interpretation of the backreaction 
effects and the form of the average equations to constraint the mean anisotropy of the accelerating 
region. To proceed we introduce the following domain-dependent, dimensionless set: 

n '-<ep ^"W Q ~~W? { ' 

and 

3A 

°» = W <17) 

We note that £l p , 0, a and £l\ (assuming A > 0) are positive definite, while and fig can 
take negative values. In particular, Qq is negative when Q is positive, which happens when the 
mean-square fluctuation of the shear is small compared to that of the expansion (see (16d)). In 
terms of these parameters, and after setting A = 0, condition (7) reads 

n p + n a + n n + n Q = 1 , (is) 

while (15) becomes 

n Q <-ln p -Q a . (19) 

Combining this result with constraint (18) we can recast the condition for accelerated expansion, 
without any Q, a and fig contributions, into the following more familiar form [8, 9] 

n„ + ffi*>§. (20) 

We may constrain the shear anisotropy of the accelerating averaged domain by combining con- 
ditions (18) and (20). For instance, setting flu = and substituting (20) into (18) gives 

o <n a < -i-n e . (21) 

This is a self-consistent condition as long as fig < —1/3, or equivalently Q > 2(0) 2 /9 (see 
definition (16d)). Following (9), the latter requires the mean-square fluctuations in the volume 
expansion to exceed those in the shear by at least 2(0) 2 /9. Note that, when Qq — > — 1/3~, 
accelerated expansion occurs when there is minimal shear anisotropy (i.e. for Q a — > 0). 



purposes. If adopted, however, the averaged Raychaudhuri equation (e.g. see (13)) will take the form 

3— = ~^<P> + § (<e 2 > - <9> 2 ) - 2S 2 - 2{Z ab * ab ) - 2(a) 2 + A , (14) 
(id ^ ' 

with an analogous expansion for (7). Clearly, when the 'background' is isotropic all the S-terms vanish identically. 
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5 Discussion 



The averaging problem in general relativity and cosmology is an issue of major importance that 
remains largely unresolved. As yet and despite the efforts there is no general consensus on what 
the approach to spatial averaging should be. Generally, averaged inhomogeneous spacetimes 
are known to evolve differently from the standard cosmological model, represented by the FRW 
universe. It is also known that spatial averaging leads to effective source terms in the Einstein 
equations, analogous to Reynolds stress terms in turbulent fluid flows, which carry what are 
commonly referred to as backreaction effects. In fact, these effective dynamical sources are 
what an observer would expect to measure. However, although the measurements take place in 
the real inhomogeneous world, their interpretation depends on the chosen fitting cosmological 
model. So far, the latter has been the Friedmann universe and the studies have shown that 
spatial averaging can seriously modify its standard dynamics, particularly its expansion rate. 
In view of the uncertainty about the averaging procedure and its physical consequences, it is 
important to interpret spatial averaging within a range wider than the FRW models. 

In this study we assumed a Bianchi-type cosmology with pressure-free matter and zero or 
isotropic 3-curvature. Our intention was to calculate the effective evolution of the averaged 
counterparts of these simple anisotropic models. Allowing for explicit average-shear terms in 
the equation, we recast the standard formulae in a closed Bianchi-type form. Thus, in our case, 
the backreaction source term is given by the difference between the mean-square fluctuations in 
the expansion and the shear scalars. For a vacuum fitting model with Euclidean 3-spaces, we 
found that the source term vanishes and the averaged equations retain their local form. When 
pressureless matter is introduced, however, the backreaction effects can dominate the early 
evolution of averaged Bianchi / models and prevent the familiar Kasner-type singularity from 
forming. Assuming axial symmetry, we have considered the 'bounce conditions' for different 
types and degrees of anisotropy. The possibility of avoiding the Kasner regime in averaged 
inhomogeneous Bianchi / universes, raises the question as to whether the same backreaction 
effects could also alleviate the Mixmaster behaviour in averaged Bianchi IX spacetimes. By 
giving a Bianchi-type form to the average equations, we were allowed to consider the limits 
that a period of backreaction triggered accelerated expansion might put on the anisotropy of 
the averaged domain. Our results show that these constraints are largely compatible with the 
observed high isotropy of the universe. 

When averaging relativistic spacetimes one should always keep in mind that there is no 
unique approach so far (e.g. see [24]-[26]). In addition to that, there is another central concern 
that besets all averaging methods and creates uncertainty about their conclusions (e.g. see [13]). 
Averaging takes an ensemble of spacetimes (Mj, gi) - all of which satisfy the standard Einstein 
equations - and introduces a procedure by which we can define an 'average' spacetime (M, g) 
associated with them. This process necessarily loses information about the (Mj, gi) which could 
lead to spurious astrophysical conclusions about the behaviour of (M, g). For instance, if every 
member of the ensemble (Mj, g{) is a dust-filled FRW universe with an initial singularity, then 
it seems very peculiar if an 'average' of these behaviours, every one of which is singular, could 
give rise to a non-singular 'average' universe. Yet, this is in effect what is happening in many of 
the studies of the averaged behaviour of FRW universes. Each member of the ensemble obeys 
the strong-energy condition and possesses decelerating expansion, but the averaged universe 
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accelerates and in a sense it effectively violates the strong energy condition. This is exactly 
the situation that is being appealed to in a number of recent studies, and also emerges in our 
examples. It remains to be seen whether the extension of the averaging procedure to deal with 
non-scalar quantities can produce a consistent resolution of this information-loss problem by 
controlling the type and the quantity of the lost information. We know that the process of 
perturbing or averaging the Einstein equations is not a sequence of operations that commutes. 
The definition of an averaged spacetime would be the result of solving the Einstein equations 
in general and then averaging the solution. In practice, we average them and then solve the 
averaged equations. This suggests that we might narrow the gap between the two results of these 
approaches by averaging an exact inhomogeneous cosmological solution of Einstein's equations 
and comparing the result with the averaging process applied to a FRW ensemble. Our analysis 
has been made in this same spirit of testing the self-consistency of an averaging procedure in 
new situations to determine whether we can rely upon the more familiar ones. 
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